Optical pumping into many-body entanglement 
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We propose a scheme of optical pumping by which a system of atoms coupled to harmonic oscil- 
lators is driven to an entangled steady state through the atomic spontaneous emission. It is shown 
that the optical pumping can be tailored so that the many-body atomic state asymptotically reaches 
an arbitrary stabilizer state regardless of the initial state. The proposed scheme can be suited to 
various physical systems. In particular, the ion-trap realization is well within current technology. 



Generating entanglement is one of the landmarks of 
modern physics, which has provoked a huge body of 
work. Entanglement — quantum correlation — is distin- 
guished from the classical counterpart by its coherent 
nature [l|. While the correlation itself can be naturally 
built up by letting multiple bodies mutually interact, 
the difficulty usually arises in retaining the coherence, 
since opening an interaction channel inevitably brings in 
adversarial environmental effects as well. Consequently, 
both the theoretical and experimental efforts have been 
devoted mostly to finding a way to isolate the system, 
whereby the problem reduces to identifying proper time 
points at which the unitary dynamics gives a useful form 
of entanglement [2| . This requires preparation of definite 
initial states and precise control of the system parameters 
and timing. In addition, the entanglement generation can 
be aided by measurements Q. 

A conceptually different approach is to tailor an open 
quantum system in such a way that its irreversible dy- 
namics drives the system into a steady state that is en- 
tangled • This approach has apparent advantages 
over those relying on the unitary evolutions. Firstly, the 
generated entangled state is robust against decoherence 
because the intrinsic decoherence channel of the system 
itself is exploited as a resource. Secondly, the steady 
state is determined by the relationship (e.g., ratio) be- 
tween different parameters, not by their exact values or 
temporal profiles, and moreover is independent of the 
initial state. This largely alleviates the otherwise nec- 
essary manipulation and control mentioned above and 
also provides extra robustness against fluctuations of the 
parameters (e.g., of laser fields). In some sense, such a 
possibility is recognizable from the fact that the envi- 
ronmental effects in ordinary quantum optical systems 
are well accounted for in terms of the zero-temperature 
bath (lo| . This means that such an environment always 
plays the role of an entropy sink that renders the system 
perfectly coherent. For single atoms, optical pumping — 
more specifically, coherent population trapping — indeed 
uses this environmental decay, i.e., the spontaneous emis- 
sion, to prepare an arbitrary pure superposition of two 
ground states [Tlj]. Generalizing this idea to multiple 



atoms, however, seems daunting and somewhat counter- 
intuitive at first sight because the spontaneous emission 
occurs independently in each atom, destroying both the 
correlation and coherence between the atoms. 

In this paper, we show that the concept of optical 
pumping can indeed be applied to the cases of many 
atoms when they are coupled to discrete quantum sys- 
tems, exemplified here by harmonic oscillators. To be 
specific, we show that as far as the geometry permits, an 
arbitrary stabilizer state (a broad class of many-body 
entangled states) of the atoms can be generated as a 
steady state through the atomic spontaneous emission. 
This contrasts with all previous works, mostly with two 
atoms, which have rather been case studies simply re- 
vealing that competition between a continuous driving 
and the environmental decay can lead to a particular en- 
tangled (but far from maximally entangled) steady state 
[Hill or have shown that atoms driven by entangled fields 
can have a steady-state entanglement [y, |7j . Our system- 
atic study, on the other hand, suggests a clear-cut mech- 
anism to generate a broad class of entangled states with a 
high fidelity using classical pumping fields. Recent stud- 
ies offer formalisms for engineering system-environment 
couplings to obtain a desired many-body entanglement 
as a steady state [l, However, as the resulting Lind- 
blad operators are multipartite ones, which require the 
interaction between the system and the bath to be ef- 
fectively many-body, the formalism does not necessarily 
make its physical realization apparent, except for some 
special instances leading to only a restricted set of en- 
tangled states [8| or those aided in large part by uni- 
tary manipulations Our scheme, on the other hand, 
deals with realistic physical models: we show that our 
idea is directly applicable to ion-trap systems and also to 
many other systems by allowing only single-atom unitary 
transformations in addition. Thanks to these exclusive 
features, our scheme would find straightforward applica- 
tions, e.g., in establishing a robust entanglement channel 
between computational nodes or continuously stabilizing 
an entangled state generated by a different means. As 
highlighted in recent literatures, our scheme would also 
have implications in the context of quantum many-body 
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FIG. 1: Involved atomic levels and transitions. g L ' R denotes 
the coupling rate between the atom and the harmonic oscil- 
lator, ui L ' R the Rabi frequency of the pumping field, 5 L ' R the 
detuning, and T L ' R the spontaneous emission rate. 



simulation and quantum computation [1, ID, [l2 , 13 1 . 

Stabilizer states represent a large set of entangled 
states including error-correcting codes and graph states 
3) TBI . Let us employ the following convention. Con- 
sider an 7V-qubit system with a 2 Ar -dimensional Hilbert 
space T-L and M stabilizers = (2>,- =1 u m with 1 < /i < 



M, where u^j € {o^ ,<r] ,a?,Ij} is one of the Pauli op- 



erators or the identity operator acting on the j-th. qubit. 
As [Sp,S„] = for every pair and (S^) 2 = 
the Hilbert space can be divided into 2 AI subspaces, 
each with dimension 2 N ~ M , according to the eigenval- 
ues of the stabilizers: H = ® Sm=±1 "H(si, sm), 
where S M |*) = |*) for |*) E H{s 1 , s M )■ Our 
aim is to pump all the population into the subspace 
H(+l,+l,..., +1). If M=N, the dimension of this Hilbert 
space is 1, i.e., it represents a single pure state, but we 
can also consider more general cases where M < N, in 
which case the steady state is a certain mixed state in an 
encoded subspace rather than a pure state. 

The stabilizer states that can be prepared in our 
scheme depend on the geometry of the arrangement of 
harmonic oscillators and atoms. For simplicity, let us 
consider a case where one harmonic oscillator is coupled 
to the atoms. Generalizing this setting to other cases 
with more harmonic oscillators is straightforward. The 
underlying idea is to have the atoms collectively coupled 
to harmonic oscillators, thereby splitting the otherwise 
degenerate energy levels in the excited manifold of the 
system according to, in our case, the parity. We then 
resolve the energy splittings so as to optically pump the 
population with a particular parity into the disjoint sub- 
space. As will become apparent, such an optical pumping 
is suitable for the generation of stabilizer states. 

We use three-level atoms with two ground levels \L) 
and \R) and an excited level \E), as shown in FIG. [U 
The two atomic transitions \L) «-> \E) and \R) o \E) 
are coupled respectively to orthogonal polarizations of 
light. Let us denote by a (a 1 ) the annihilation (creation) 
operator for the harmonic oscillator. We assume that the 



harmonic oscillator interacts with the atom in such a way 
that they exchange energy, which can be described in the 
interaction picture by the following Hamiltonian: 

N 

ffah - £ a (gf \E) j (L\ + gf \E) j (R\) + H.c, (1) 
3=1 

where the subscript j represents the j-th atom. The cou- 
pling rates gf and g R are assumed to be independently 
adjustable. On top of this coupling, we also apply classi- 
cal pumping fields with Rabi frequencies ujf and lo r and 
detunings S L and S R , respectively, as in FIG.[TJ For con- 
venience, let us call the former an L-ficld and the latter 
an i?-field. This driving can be described in the interac- 
tion picture by the following Hamiltonian: 



H n = 



N 
3=1 



ujfe* d 'I^WLl+wfe* 4 \E) . (R\ + H.c. (2) 



Once the atom is excited, spontaneous emission to \L) 
and \R) takes place with rates T L and T R , respectively, 
which can be described by the following master equation: 



N 



dr 



ah 



j=l x=L,R 



(3) 



where p is the density matrix of the system and £J (p) = 
T-(\x) j (E\ p \E) 3 (x\ - I \E). (E\ p - \p \E), (E\) . 

To begin with, suppose only two atoms are coupled to 
the harmonic oscillator and only the transition \R) -H> \E) 
is coupled, i.e., g\ = g% = and g R = g R = g. We can 
denote the state by \(j>, n), where \<j>) denotes the atomic 
state and |n) the n-th excited state of the harmonic oscil- 
lator. Note that the coupling between the atom and the 
harmonic oscillator does not change the number of to- 
tal excitations n t = a^a + J^j i.e., [if a h,n t ] = 0. 
Putting the pumping fields and the spontaneous emission 
aside, the Hilbert space is thus divided into mutually de- 
coupled subspaces according to n t . Furthermore, due to 
the coupling, the degeneracy in each excited manifold is 
lifted. In the manifold of n t = 1, transitions \LR, 1) «-» 
\LE,0) and \RL,1) O \EL,0) occur with coupling rate 
g, whereas transition \RR,1) O A=(\ER,Q) + \RE,0)) 

occurs with coupling rate y/2g as the atoms collectively 
interact with the harmonic oscillator. Consequently, as 
shown in FIG. [21 a), the excited manifold shows energy 
splittings, whose amounts are determined by the number 
of atoms in states \R) and \E). Formally, one can define 
n s = J2j \R)j (R\ + \E)j (E\ and then the energy splitting 
is given by ±y/n^g (except for uncoupled states). 

We resolve these energy splittings by adjusting the de- 
tunings to achieve selective pumping of the population. 
As we aim to generate an entangled state in the ground- 
state manifold (i.e., n t = 0) and the spontaneous emis- 
sion always projects the population to the lower man- 
ifold, we need to consider only the two lowest mani- 
folds of nt = and n t — 1. FIG. [^b) depicts the 
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FIG. 2: (a) Energy splitting in the manifold of one excitation 
for two atoms (normalization is ignored), (b) Optical pump- 
ing from the odd-parity subspace to the even-parity subspace. 
Straight arrows represent classical pumping fields and wiggly 
arrows spontaneous emissions. ui r /\/2 and u /2 are the cor- 
responding effective Rabi frequencies. 



pumping process from the subspace with of erf = — 1 
to that with erf erf = +1. Here, the Pauli operators 
are defined over the ground levels of the atoms, e.g., 
a? = \L) j (L| - \R) j (R\. Note that the L- and i?-fields 
couple the two manifolds in different ways. While the 
L-field induces transition in such a way that when the 
state is excited n s is increased by one, the i?-field pre- 
serves n s . The spontaneous emission, on the other hand, 
only decreases th by one while either decreasing n s by 
one or preserving n s (see FIG. HJb)). The overall effect 
of the pumping is thus that the L-field increases n s by 
one while the i?-field decreases n s by one. Thanks to the 
energy splittings, one can apply these optical pumpings 
selectively to the population with a particular n s : in or- 
der to pump the population with n s to the subspace with 
?i s + l (n s — 1), we apply an L-field (i?-field) to every atom 
with detuning ±^/n s + lg (zt^/n^g). 

It is easily seen that the above pumping process drives 
a two-qubit system into the subspace with <S M = erf erf = 
+1. Generalizing this idea to the case of an arbitrary 5^ 
is straightforward. Suppose first that we need to pump 
the population with of of • • • ojy = —1 to the subspace 
with cf erf • • -erf = +1. We achieve this by pumping 
every population with an odd n s , i.e., that with an odd 
(N + erf + • • • + Ojy)/2, to the subspace with n B ± 1. For 
this, we apply to every atom L-fields with detunings V2g, 
\fig, V&g, and so forth. In the same fashion, we apply 



to every atom i?-fields with detunings —g, —\/3g, —y/5g, 
and so forth. Here, we exploit red (blue) detunings for In- 
fields (R- fields) , but this is not mandatory. For a general 
5^, we can achieve the same kind of optical pumping 
by rotating the basis of the operation for each j-th atom 
according to u^j . If u^j = Ij , we turn off all the couplings 
and pumping fields for the j-th atom so that the atom 
is excluded from the pumping process. Otherwise, we 
adjust the coupling rates so that — a*g and g^ — 
(3*g, where ctj \L) ■ + /3j \R)j is the +1 eigenstate of m w -, 
and instead of an 7?-field with Rabi frequency lu^, we 
apply an L-field with a^uf and an R- field with 
[16j| . Likewise, instead of an L-field with Rabi frequency 
Wj, we apply an L-field with PjUjj' and an R- field with 

A general stabilizer state is now generated by per- 
forming the optical pumping processes for all S^'s one 
by one as described above. An important question here 
is how the performance, namely, the characteristic time 
for the system to reach the steady state, scales as the 
number of atoms increases. Unfortunately, it turns out 
that the time increases exponentially with the num- 
ber of stabilizers M. In order to see this, it is con- 
venient to think of the pumping process as a random 
walk (or hopping) among the 2 M subspaces H(si, sm)- 
For example, suppose a linear cluster state \^>n) — 

2- N ' 2 ®f- 1 \\L) J (L\ + \R) J (R\af +1 )®ti(\L) k + \R) k ) 
of four atoms (N = 4), stabilized by <Si = erf- erf I3I4, 
S 2 = ofa?cr£h, S 3 = Iicf erf of , and S 4 = hl 2 a$a$, 
and consider the pumping step for S\ . The original aim of 
this step is to pump the population in H{— 1, S2, S3, S4) 
to H(+l , S2, S3, S4). However, as it destroys the coher- 
ence of atom I and atom 2, the values of 1S2 and £3, 
which contain a Pauli operator on atom 1 or atom 2, are 
also affected. Consequently, the pumping for Si in fact 
occurs from %{— 1, s%, S3, S4) to one of H(+l, ±1, ±1, S4) 
randomly. The performance is then assessed by estimat- 
ing the average number of steps to reach the subspace 
H(+l , +1) by such a random walk, along with the ac- 
tual pumping time for each step. The latter increases 
as the number of Pauli operators in iS M , which we de- 
note by k^, increases because the energy gap to resolve 
g{^fk^ — y/k^ — 1) decreases. This scales only polyno- 
mially and can be bounded as long finite for 

all S^, which is usually the case. However, the former, 
the number of steps to reach the steady state, increases 
exponentially as < M(2 /m ) M , where each pumping step 
affects at most l m other stabilizers. Our scheme is thus 
more relevant for a moderate number of atoms. 

FIG. [3] shows the results of our simulation for the gen- 
eration of linear cluster states, which is done by numeri- 
cally integrating the master equation ([3]). As the pump- 
ing exploits the energy splitting, both the pumping field 
and the spontaneous emission should be weak enough 
compared to the energy gap so that the energy split- 
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FIG. 3: Optical pumping to linear cluster states for N = 
2, 3, 4 when a fully mixed state is taken as the initial state for 
(a) A = 0.05 and (b) A = 0.02. A is a constant determining 
the degree to which the energy splitting is resolved (see text). 



ter each pumping and detune the transition by an ac 
Stark shift for those corresponding to Ij . It is easily seen 
that this modified scheme is well suited to systems of 
atoms sitting in a common cavity mode. Another modi- 
fication can be made by noting that the harmonic oscil- 
lator can in fact be replaced by a two-level atom, where 
the interaction between atoms is the spin-exchange in- 
teraction. Along with single-qubit manipulations as de- 
scribed above, which are relatively easy to do, our scheme 
can thus be applied to various Heisenberg spin systems 
realized with atoms, where it is understood that auxil- 
iary levels are needed to represent another ground level 
\L) and to realize spontaneous emissions effectively if re- 
quired 21, 22l |. 



tings can be resolved, 
ters as r = \g(^fk~l - 



For this, we choose the parame- 

r L = r 



i 



T/2, and 



lj = w = r / where A < 1 is a small constant, 
and each round of the pumping for one stabilizer is per- 
formed for a period of time tt/T. Our particular choice of 
the parameters is, however, not mandatory as the pump- 
ing is by nature insensitive to the exact values of the 
parameters. Taking a fully mixed state of atoms along 
with the ground state of the harmonic oscillator as an 
initial state, FIG. [3] shows the time evolution of the fi- 
delity of the state to the desired linear cluster state. The 
weaker the pumping field and the spontaneous emission 
are (smaller A), the higher is the final fidelity, but the 
longer it takes to reach the steady state. One could vary 
the parameters in time to optimize the performance, e.g., 
by taking larger A during the initial transient period and 
taking smaller A afterward to obtain a higher fidelity. 

It is evident that the present scheme is well suited to 
the existing ion-trap systems used for quantum informa- 

For example, Ca + ions can be 
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tion processing 
used, where the two hyperfine levels of S1/2 represent 
\L) and \R), respectively, and one of D5/2 represents \E) 
[l9(. The coupling between each ion and the center of 
mass mode of the ions g^ ' can be controlled individu- 
ally by focused beams. As D 5 / 2 is metastable, the spon- 
taneous emission rate T L,R can be controlled by adjusting 
pumping fields via P3/2 (and repumping fields via P\/2 to 
make the transition closed) . As the typical coupling rate 
is 9 ^ 100 kHz, the system reaches the steady state in 
a few hundred milliseconds for a small number of atoms, 
as can be inferred from FIG. [3] One can also incorpo- 
rate sympathetic cooling into the pumping sequence to 
overcome motional heating [20(. 

The present scheme with slight modifications can be 
applied to many other physical systems. First, we can 
consider a case where only one of the two atomic tran- 
sitions is persistently coupled, i.e., g L = and g R 7^ 0. 
This enables the optical pumping for stabilizer ■ a? . 
For different stabilizers, one performs single-qubit rota- 
tions to the atoms corresponding to a ■ : before and af- 
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